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We study the light scattered from randomly rough, one-dimensional self-affine fractal silver sur- 
faces with nanoscale lower cutoff, illuminated by s- or p-polarized Gaussian beams a few microns 
wide. By means of rigorous numerical calculations based on the Green theorem integral equation 
formulation, we obtain both the far- and near-field scattered intensities. The influence of diminishing 
the fractal lower scale cutoff (from below a hundred, down to a few nanometers) is analyzed in the 
case of both single realizations and ensemble average magnitudes. For s polarization, variations are 
small in the far field, being only significant in the higher spatial frequency components of evanes- 
cent character in the near field. In the case of p polarization, however, the nanoscale cutoff has 
remarkable effects stemming from the roughness-induced excitation of surface-plasmon polaritons. 
In the far field, the effect is noticed both in the speckle pattern variation and in the decrease of the 
total refiected energy upon ensemble averaging, due to increased absorption. In the near field, more 
efficient excitation of localized optical modes is achieved with smaller cutoff, which in turn leads to 
huge surface electric field enhancements. 

PACS numbers: 290.5880, 290.4210, 240.6680, 180.5810. 



I. INTRODUCTION 

Since the early days of fractality,^] the scattering of 
electromagnetic (EM) waves from fractal surfaces has 
been a field of intense activity. Physical fractals ap- 
pear ubiquitously in nature, possessing fractal properties 
within a broad, however finite, range of scales. Therefore 
the study of classical wave scattering from fractals is a 
problem of interest not only from a fundamental point of 
view, but also from the practical knowledge that probing 
technologies such as surface optical characterization, re- 
mote sensing, radar, sonar, etc. can yield about a wide 
variety of systems. In fact, it is now well understood 
that many naturally occurring surfaces exhibit scale in- 
variance, particularly in the form of self-affinity. H, |^ 

There exists a large amount of theoretical works de- 
voted to wave scattering from fractal surfaces. 0, | ^ , |^, |^, 
|, H 0, 0, g, [16[ 0, H, |0[ iJlf Most 
of them make use of approximations, such as Kirchhoff 
approximation (KA) or perturbation methods, in order 
to obtain analytical expressions that are useful in cer- 
tain regimes, but thereby imposing a constraint on the 
length scales over which fractality might be present. Oth- 
ers include calculations within the Rayleigh hypothesis 
and/or restricted to perfectly conducting surfaces. And 
only very recent works are capable of dealing with arbi- 
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trarily rough metal or dielectric surfaces, fri 
on the basis on the Green's theorem integral equation 
formulation for rough surface scattering. ||2^, |2j, |2J] 

On the other hand, in addition to far field scattering, 
the near EM field on self-affine fractals has attracted 
a great deal of attention. ||, |7[ |2|, |2|, |l], |3|] 
To a large extent, the interest lies in the observation 
in the near field intensity distributions, through photon 
scanning tunneling microscopy (PSTM), of large con- 
centrations of EM field intensity on bright spots) . l3l|, 
|3^ , HJ, m ^ The appearance of such local- 

ized optical modes is mediated by the roughness-induced 
excitation of surface-plasmon polaritons (SPP) on very 
rough surfaces with subwavelength features. The oc- 
currence of large surface EM fields is indeed crucial to 
the EM mechanism in Surface-Enhanced Raman Scat- 
tering (SERS) and to other non-linear surface optical 

Furthermore, the 



processes. [E6l 



30 31 



predicted, enormously high values of the surface electric 
field intensity have been related to the recent observation 
of SERS single molecule probing. |4|, H 

It is our aim to employ the above mentioned rigor- 
ous Green's theorem formulation to study the far-field 
and near-field scattered from one-dimensionally rough, 
self-affine fractal Ag surfaces, restricted to one dimension 
for the sake of computational limitations. The scatter- 
ing model is described in Sec. II. We focus on fractals 
preserving self-affinity from tens of microns to subwave- 
length dimensions. Actually, we analyze the influence of 
the lower scale cutoff as is reduced from 100 nm (as 
used in Refs. |lj, 18, ^ ^ ^2|) to a few nanometers. It 
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should be mentioned that the effect on the far field of the 
lower scale cutoff in deterministic Koch fractals has been 
previously studied. The results for the far-field and 
near-field calculations are presented in Sees. Ill and 
respectively. Finally, the main conclusions derived from 
this work are summarized in Sec. 0. 



II. SCATTERING MODEL 

The scattering geometry is depicted in Fig. |^. A 
monochromatic, linearly polarized Gaussian beam of fre- 
quency Lo and 1/e amplitude half width W impinges at 
an angle on a randomly rough surface z — (^{x). The 
rough interface separates vacuum from a semi-infinite sil- 
ver volume occupying the lower half-space [z < C{x)] and 
characterized by an isotropic, frequency-dependent ho- 
mogeneous dielectric function e{uj). 

With the aim of solving this scattering problem for 
arbitrarily large roughness parameters, we make use of 
the scattering integral equations based on the the appli- 
cation of Green's second integral theorem, [p^ It is well 
known that this full vectorial formulation is considerably 
simplified when restricted to ID surfaces and linear po- 
larization, thereby being reduced to four scalar integral 
equations for the only nonzero component of either the 
electric field amplitude 



with 



H^") EE Hi'^ (r, w)x + iJi") (r, w)z, 
for s polarization, or the magnetic field amplitude 
H(p) = H^P\r,uj)y, 

with 

E(p) EE E^P^ (r, w)x + EiP'^ (r, w)z, 

for p polarization. H^] Two of these surface integral 
equations can be used as extended boundary conditions, 
which (recalling the continuity relations across the inter- 
face) lead to a system of two coupled integral equations 
for the surface field and its normal derivative. This can 
be solved numerically in the form of linear equations by 
means of a quadrature scheme. |2^, ^ Once these source 
functions are obtained, it is straightforward to calculate 
the far field scattered intensity, |23] and in general the 
electric (magnetic) field amplitude for s (p) polarization 
at any point in vacuum or inside the metal from two of 
the integral equations. ||29| In addition, from the latter 
integral equations, the corresponding expressions for the 
magnetic (electric) field amplitudes for s (p) polarization 
can be simply obtained by making use of Maxwell equa- 
tions (cf . Ref . ^ in the case of p polarization) . These 
expressions can be useful in near electric or magnetic 



field calculations. However, they exhibit non-integrable 
singularities when approaching the surface profile. 

In this regard, the calculation of the surface magnetic 
field intensity for s polarization has been done in a simple 
manner by exploiting the connection between its normal 
and tangential components on the vacuum side 
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with the surface electric field and its normal derivative 
from the vacuum side. 
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with 7 = [1 + [C'(a;)]^]^/^. The corresponding expressions 
for the surface electric field components for p polarization 
are given by: [E2| 
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As mentioned above, E^'\x), F^'l (x), H^P^x), L^p^x) 
constitute the source functions of the integral equations 
in this scattering configuration. 

As a model describing many naturally occurring sur- 
face growth phenomena exhibiting self-affine fractality, 
we have chosen that given by the trace of a fractional 
Brownian motion through Voss' algorithm, ||, ^ In 
addition, this model yields self-affine fractal structures 
that resemble fairly well the properties of some SERS 
metal substrates. Q The ensembles of realiza- 

tions thus generated are characterized by their fractal 
dimension D = 2 — 71 [Ti being the Hurst exponent) and 
rms height 5. (In order to avoid the inherent ambiguity 
in the definition of the rms height for self-affine fractals, 
5 refers in our calculations to the rms height defined over 
the entire fractal profile with length Lf — 51.4 /im. Re- 
call that 5 depends on the length Ax over which it is mea- 
sured through 5 — l^~^Ax^, I being the topothesy. 
Thus the topothesy I will be also given. From each gen- 
erated fractal profile with Nf points and length Lf, se- 
quences of N points (with constant N /Nf) are extracted 
to obtain similar profiles with identical properties except 
for the lower scale cutoff as determined by — Lf/Nf. 
(Strictly speaking, should be given by the minimum 
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length scale above which the ensemble of such realiza- 
tions exhibit self-afSnity, resulting in a value typically 
larger than the mere discretization cutoff. ^ ) 

From now on, we will focus below on the fractal dimen- 
sion D = 1.9. This choice is justified for we expect the 
influence of decreasing to be more significant for larger 
fractal dimensions. And recall that the effect of varying 
D has been already studied on the far field jlSj and the 
near field albeit for a relatively large lower cutoff 
^L- The values of the lower scale cutoffs hitherto consid- 
ered are = 51.4,25.7,12.85, and 6.425 nm, resulting 
from sequences of iV = 102,205,410, and 819 points ex- 
tracted from profiles with Nj = 1024,2048,4096, and 
8192 points. The length of all realizations is thus L = 
Lf{N/Nf) — 5.14/im. Actually, the final number of sam- 
pling points per realization used in the numerical calcu- 
lations is significantly higher for the sake of accuracy: 
Np = TiiN as obtained by introducing Ui =4-10 cubic- 
splined interpolating points, the latter being chosen on 
the basis of numerical convergence tests. 



III. FAR FIELD 

First, we present in Fig. ^the angular distributions of 
the speckle pattern intensities scattered from self-affine 
Ag surface profiles with identical properties (D = 1.9, 
L = 5.14 /im, and 5) except for = 51.4, 25.7, 12.85, and 
6.425 nm, as mentioned in the preceding section. Both 
s- and p-polarized incident {9q — —10°) Gaussian beams 
are considered with wavelength A = 27rc/a; = 629.9 nm 
and W = (L/4)cos6'o, the dielectric constant of Ag at 
this frequency being[^ e = —15-1-1. The upper plots 
correspond to moderately rough surfaces with 5 = 51.4 
nm and topothesy I = 24 nm (the corresponding surface 
profiles will be shown in Fig. ^), beyond the reach of 
standard approximate treatments such as the KA or per- 
turbation theories. No significant differences are found 
(note the semi-log scale in order to enhance them) for 
s polarization, nor even away from the bright specular 
peak. In the case of p polarization, a weak decrease in 
the specular peak with decreasing is barely observed, 
whereas small differences appear in the speckle pattern 
away from the specular peak. 

We now consider rougher surfaces with 6 = 257 nm 
and topothesy I — 143 nm (see Figs. ||c and ||d). Ap- 
preciable changes are observed in the s-polarized speckle 
with decreasing ^l, although the patterns are qualita- 
tively similar. However, in the case of p-polarized speckle 
patterns, differences are both quantitative and qualita- 
tively notorious. Note the absence of specular peaks in 
both polarizations. Therefore, in light of Figs. ^ and 
it is seen that the details at nanometer scales produce no 
change in the s-polarized speckle patterns for moderately 
rough self-affine profiles, appearing only small quantita- 
tive variations for very large rms height. For p polariza- 
tion, nonetheless, small changes are already observed for 
moderately rough profiles, becoming dramatic for large 



roughness. 

Next, the influence of nanoscales on the mean scat- 
tered intensities is analyzed. This is done in Figs. || after 
averaging speckle pattern distributions over an ensemble 
of Nrea — 200 realizations. Experimentally, this can be 
achieved either by also ensemble averaging, or by illumi- 
nating a very large surface and collecting a large amount 
of speckles at every scattering angle (or a combination of 
both). In the case of the self-afhne fractals with 6 = 51.4 
nm (see Figs. ^ and ^), the diffuse component of the 
mean scattered intensity in s polarization is not altered 
at all with decreasing lower scale cutoff. As for the specu- 
lar, which amounts to 32%, and not shown in Fig. |^), 
the same is true. The p-polarized diffuse component in 
Fig. IIId shows no significant variation at angles away from 
the specular direction, whereas a very small, but ob- 
servable, decrease in the specular direction is observed 
when the lower scale cutoff diminishes. On the other 
hand, the p-polarized specular component, not shown in 
Fig. ||b, exhibits a larger decrement (from ~ 32% down 
to ^ 24%) with diminishing ^l. Interestingly, despite the 
specular components are not included in Figs. |^a and ^d, 
the resulting angular distributions of diffusely scattered 
intensities present pronounced peaks at the specular di- 
rections. These diffuse specular peaks have been already 
discussed in Ref. |l^ for self-afhne fractals with various 
fractal dimensions but a higher lower scale cutoff; the 
results in Figs. ||a and demonstrate that such peaks 
are not substantially affected by the fractal scales below 
a hundred nanometers. 

Ensemble averaging of the speckle pattern scattered 
from very rough (S = 257 nm) self-affine fractal sur- 
faces also tends to wash out differences arising from 
nanoscale fractal details (see Figs. ^ and ||d). In fact, 
no differences are found in the angular distribution of s- 
polarized, mean scattered intensities (leaving aside some 
spurious speckle noise due to the limited number of re- 
alizations employed), even though quantitative changes 
have been found above in the speckle patterns for single 
realizations (see Fig. ^). Furthermore, the remarkable 
qualitative and quantitative changes in the p-polarized 
speckle spots seen in Fig. ||d reduce to significant, though 
only quantitative, variations upon ensemble averaging as 
shown in Fig. And if the effect of absorption, defined 
'AS A — 1 — S {S being the normalized reflectance) is 
taken into account, the resulting renormalized angular 
distributions of p-polarized mean scattered intensities, 
{I{9))/S, exhibit little differences with decreasing frac- 
tal lower scale cutoff (see Fig. ^) . Incidentally, note that 
the remaining speckle noise in Figs. |^ and ^ hinders the 
possible observation of a backscattering peak due to the 
multiple scattering of SPP. ||l^ 

IV. NEAR FIELD 

We now turn to the investigation of the influence on the 
near EM field of the lower scale cutoffs of the self-afhne 
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fractal surfaces whose far-field scattering properties have 
been studied in the preceding section. Near field inten- 
sity distributions are relevant for the information they 
provide on the scattering process, and can be measured 
through near-field optical microscopy. 

First, we calculate the intensities of all the EM field 
components on the surface. The results for the mod- 
erately rough fractals with S = 51.4 nm employed in 
the speckle pattern calculations in Fig. || are shown in 
Figs. |5| and Ig for s and p polarization, respectively. The 
corresponding surface profiles are shown in the bottom, 
Figs. ^ and Figs. ||d. The nonzero components of the 
EM field being plotted are (cf. Sec. II): the tangential, 
perpendicular to the incident plane, electric (respectively, 
magnetic) field and the normal and tangential (in the 
plane of incidence) magnetic (respectively, electric) field, 
in the case of s (respectively, p) polarization. For the 
sake of clarity, only the central part of the illuminated 
surface is shown. 

Before analyzing the spatial distributions, let us re- 
call what is expected within the KA, namely, the sum of 
the incident and the specularly (locally) reflected fields. 
Since the angular variation of the Fresnel coefficients for 
metals at this frequency is small, the local variations of 
the specular field can be neglected. Therefore, the KA 
field intensities are approximately given by those for a 
planar surface, namely: 



^^{s,KA) |2^ 
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for s polarization; and 
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for p polarization, TZg and TZp being the corresponding 
Fresnel coefhcients for Oq on a planar metal surface. The 
latter field intensities are plotted in Figs. ^ and ||. As 
expected, the KA does not hold even for the moderately 
rough surface profiles used therein, nor does perturba- 
tion theory (recall that the planar surface field intensi- 
ties can be considered the zeroth order approximation in 
the small- amplitude perturbation expansion of the field), 
but they provide the background about which the actual 
surface EM field intensities strongly vary. 

In the case of s polarization. Fig. ^, significant vari- 
ations appear in the surface EM field upon decreasing 
the lower scale cutoff. No appreciable differences have 
been found in the corresponding far field speckle patterns 
in Fig. ^ and ||b, so that the surface field nanoscale 
fluctuations are mostly due to the evanescent compo- 
nents. Fluctuations about the background distributions, 
Eqs. (|^), become narrower and steeper the smaller is 
^L, namely, the smaller are the surface features. This 



effect is more pronounced for the normal magnetic field, 
for which the expected KA background is indeed very 
small. Interestingly, it should be noted that the elec- 
tric field intensity resembles with positive contrast the 
surface profile (the same is true for the total, tangential 
plus normal, magnetic field intensity), which is not the 
case in general. Although not valid from a quantita- 
tive standpoint, small amplitude perturbation theory can 
provide a qualitative explanation. The first-order term of 
the surface electric field in powers of the surface height, 
which in turn gives the first-order correction to the planar 
surface background, p6|, ^ can be cast for this particular 
scattering geometry, polarization and near normal inci- 
dence in the form of a convolution integral involving the 
surface profile function. We have verified, though not 
shown here, that by simply increasing the angle of inci- 
dence up to do — 40° the resemblance is slightly lost. 

Upon illuminating with p-polarized light, see Fig. |^, 
the surface EM field intensity distributions become more 
complicated, with larger variations from one realization 
to another with diminishing ^l, and no resemblance 
whatsoever with the surface profiles. These variations 
are considerably larger than those for the far-field speckle 
patterns, indicating the crucial role played by the evanes- 
cent components. Furthermore, the excitation of SPP 
propagating along the surface and reradiating into vac- 
uum mediates the scattering process. In fact, the spatial 
frequency of the large oscillations about the background 
in the surface EM field intensity is related to the SPP 
wavelength; 1 2^ this is more easily observed in the surface 
magnetic field in Fig. ^ (also in the total electric field, 
not shown here) for the profiles with higher On the 
other hand, the surface normal electric field component 
for the smaller (see Fig. ^) reveals the appearance 
of narrow and bright spots where an enhancement of the 
field intensity of nearly two orders of magnitude occur, 
despite the relatively low value of S. This clearly mani- 
fests the crucial role played by the lower scale cutoff in 
the excitation of localized optical modes. Q 

What if the surface roughness parameter is increased? 
In Figs. ^and||, the surface EM field intensities are plot- 
ted for s and p polarization, respectively, as in Figs. |^ 
and ^, but for 6 = 257 nm. Due to the large differ- 
ences introduced by ^l, a semi- log scale is used and 
only a small surface region about one wavelength long 
is shown. In the case of s polarization, the differences 
in the surface EM field intensities from one realization to 
another with lower are remarkable, and quantitatively 
larger than those differences shown above for smoother 
surface profiles; the qualitative behavior is however sim- 
ilar. Lower-^i surfaces lead to surface EM field intensi- 
ties with higher spatial frequency components, but the 
smoothly varying envelopes correspondig to smaller spa- 
tial frequencies are largely preserved. The former are 
mostly evanescent components that contribute little to 
the far-field intensities, which are basically originated in 
the propagating components, and thus quite similar for 
all surface profiles (see Fig. ^). It should be remarked 
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that no localized optical excitations seem to appear, as 
expected since no SPP can be excited, nor significant 
EM field enhancements are found. Note that, thanks to 
the semi-logarithmic scale used in Figs. for the surface 
field intensities, a weak qualitative resemblance of posi- 
tive contrast with the surface profiles can be seen. 

The p-polarized surface EM field intensities presented 
in Fig. m reveal very strong changes in all field compo- 
nents with lower These changes are considerably 
stronger than those for the smoother surfaces (see Fig. H). 
As mentioned above, narrow, large peaks tend to con- 
centrate at bright spots in the form of localized optical 
modes, exhibiting very large electric field enhancements 
preferentially polarized along the normal to the surface, 
in agreement with the SPP polarization, although sig- 
nificantly large enhancements of the tangential electric 
field intensity are also found. In addition, large magnetic 
field intensities abound, not necessarily at the same elec- 
tric field bright spots. The polarization of the modes 
has been discussed in Ref. ^ for rougher self-affine frac- 
tals with large, but yet subwavelength, S^l. It should be 
emphasized that the nanoscale lower cutoff largely con- 
tributes to the building and strengthening of such optical 
excitations. 

Let us plot the EM field in the near vicinity of the 
rough interface. This is done in Figs. || and |l^, where 
the near-field intensity maps in a logarithmic scale for 
s and p polarization, respectively, are shown in a region 
around the origin of the surface profile with = 12.85 
nm in Figs. ^ and ||d. All the EM field components 
are included for the sake of completeness; recall that it 
has been recently reported that the magnetic field inten- 
sity can be also probed through PSTM for certain ex- 
perimental configurations. |49|] The actual surface profile, 
though not explicitly depicted, can be inferred from the 
fairly black, metal regions due to the evanescent behav- 
ior of the EM fields inside metal with a small skin depth 
{d « 25 nm). 

Note that no bright spot whatsoever is observed in the 
EM field intensity maps in the case of s polarization. 
Fig. ^ which exhibits no relevant features in the near 
field. The actual interface appears slightly blurred, since 
the minima have been reduced to enhance contrast. In- 
cidentally, the continuity of the (tangential) electric field 
is satisfied in Fig. 9 well as the (dis)continuity of the 
(normal) tangential magnetic field in Figs. |c and |d. 

On the other hand, a very bright spot is found at the 
local maximum in the central part of the intensity map 
in Fig. |lo|a: The intensity field enhancement at such spot 
is \ E\^ I \ E^'^^ p~ 10'* (recall that the maximum scale of 
EM field intensities in Fig. || is about 3 orders of magni- 
tude lower than those in Fig. |o|) . Note that there seems 
to be a bright, though weaker, magnetic spot associated 
with the optical mode (see Fig. p^b), which is nonetheless 
slightly shifted to the right of the electric field maximum, 
on the surface local minimum nearby. The electric field 
intensity decays rapidly upon moving into vacuum away 
from the bright spot, faster than expected for the evanes- 



cent decay of SPP propagating on a plane. This may 
help to explain why localized optical modes experimen- 
tally observed through PSTM yield considerably smaller 
enhancement factors, leaving aside the fact that no 
direct comparison with theoretical calculations for the 
actual experimental surface profile are available. With 
regard to the electric field orientation on the bright spot. 
Figs, and [l0| d indicate that the normal electric field 
component is resTionsible for the bright spot, in agree- 
ment with Figs, gb and |^c; this component corresponds 
to I i?z P in regions near surface maxima and minima (see 



Fig. Kd) and to\E^f near vertical surface walls (see 
Fig. 10:). It can be also observed that, as expected, the 



normal electric field is discontinuous across the interface, 
the tangential component being continuous. The conti- 
nuity of the (tangential) magnetic field for p polarization 
is evident in Fig. |lo|a. 

The statistics of the surface electric field enhance- 
ment has been studied elsewhere,}^ confirming a sig- 
nificant increase with decreasing It should be re- 
marked that such large surface electric fields plays a 
decisive role in SERS and other surface non-linear op- 
tical processes; p^, p9| , |39| | in particular, local values 
of ^ 10^ can help to explain the existence of bright 
spots at which SERS single-molecule detection has been 
claimed, iol, Ell El 



V. CONCLUDING REMARKS 

To summarize, we have studied the influence of the 
nanoscale (in the region below a hundred nanometers) 
lower cutoff on the scattering of light from one- 
dimensional, self-affine fractal Ag surfaces with large 
fractal dimension D = 1.9 and for both moderately (5 = 
51.45 nm) and large {b — 257 nm) surface height devia- 
tions. Since no approximate methods are applicable for 
such roughness parameters, numerical calculations have 
been carried out on the basis of the Green's theorem in- 
tegral equation formulations, extended to account for all 
EM field components in the near field region. Both far- 
field and near-field distributions have been obtained for 
s and p polarization (no depolarization takes place in 
this scattering geometry) . We have presented far- field in- 
tensity distributions corresponding to single realizations 
(speckle patterns) and averages over an ensemble of re- 
alizations. 

It has been shown that the details corresponding 
to scales below a hundred nanometers in moderately 
rough self-affine profiles produce no change in the far 
field speckle pattern, appearing only small quantitative 
variations for very large rms height. For p polariza- 
tion, nonetheless, small changes are already observed for 
moderately rough profiles, becoming dramatic for large 
roughness. Ensemble averaging tends to wash out the 
differences. In the case of s polarization, the angular 
distributions of mean scattered intensities are practically 
identical, not only for 8 = 51.45 nm, but also for 8 = 257 
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nm. For p polarization, the differences encountered in the 
speckle patterns are largely suppressed in the mean scat- 
tered intensity distributions; the only effect of the vary- 
ing nanoscale cutoff manifests in the total reflectance 
(which is smaller for decreasing due to larger absorp- 
tion). 

Nevertheless, nanoscale features have a strong im- 
pact on the near EM field distributions due to the rele- 
vant role played by the evanescent components. The s- 
polarized surface EM field intensity exhibits oscillations 
with higher frequency the smaller is. The amplitude 
of such oscillations, and thus the changes from one pro- 
file to another with different S^l, increases with the rms 
height S. No significant electric field enhancements are 
found in this polarization. Incidentally, the surface and 
near electric and magnetic field intensities vaguely resem- 
ble the surface profile with positive contrast despite the 
large surface roughness of the profiles. In the case of p po- 
larization, drastic changes with decreasing are found 
in the surface EM field stemming from the roughness- 



induced excitation of SPP. Large, narrow peaks {local- 
ized optical modes) tend to appear with either increasing 
rms height or decreasing nanoscale, predominantly en- 
hancing the normal electric field intensity (as expected 
from the SPP electric field orientation), and rapidly de- 
caying into both vacuum and metal. Our near field in- 
tensity maps around optical excitations can help to inter- 
pret PSTM experimental results. In addition, the large 
field enhancements associated with such localized optical 
modes on self-affine fractal metal surfaces can shed light 
onto SERS and surface nonlinear optical phenomena. 
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FIG. 1: Illustration of the scattering geometry. 
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FIG. 3: Same as in Fig. gbut for the mean scattered intensity 
averaged over Nr — 200 realizations. 
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FIG. 2: Angular distribution of the speckle pattern intensity 
scattered from Ag fractal surface profiles with L — 5.14 /im, 
D = 1.9, and S — 51.4 nm (top, a-b) and S — 257 nm (bottom, 
c-d), difltering only in = 51.4,25.7,12.85, and 6.425 nm. 
illuminated by a Gaussian beam with 9o — —10°, A — 629.9 
nm, and W = L/4cos^o- Left (a and c): s polarization; right 
(b and d): p polarization. 
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FIG. 4: Same as in Figs. |^ and ^ but renormalizing by the 
total reflectance S. 
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FIG. 5: Surface electric (a) and magnetic, normal (b) and 
tangential (c), field intensities on the center spot of the illu- 
minated area L = 5.14 /xm for s-polarized scattering with Oo = 
—10°, A = 629.9 nm, and W = L/4 cos 6o, from Ag fractal sur- 
faces with D = 1.9, S = 51.4 nm, and = 51.4,25.7, 12.85, 
and 6.425 nm. The KA field intensity is also included (see 
text), (d) The corresponding surface profiles. 
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FIG. 6: Surface magnetic (a) and electric, normal (b) and 
tangential (c), field intensities on the center spot of the illu- 
minated area L = 5.14 nm for p-polarized scattering with = 
-10° , A = 629.9 nm, and W = L/4 cos 6*0, from Ag fractal sur- 
faces with D = 1.9, S = 51.4 nm, and = 51.4,25.7, 12.85, 
and 6.425 nm. The KA field intensity is also included (see 
text), (d) The corresponding surface profiles. 
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FIG. 7: Surface EM field intensities (in a semi-log scale) for 
s polarization as in Fig. ^ but for rougher surface profiles (d) 
with S — 257 nm. 
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FIG. 9: Near field intensity images in s polarization [(a) elec- 
tric, (b) magnetic, the former split into (c) x, horizontal com- 
ponent, and (d) z, vertical component, all of them in a logjQ 
scale! in an area of 0.5x0.5 fim^ close to the fractal surface in 
Fig. Md with S = 257 nm and = 12.85 nm. Other parame- 
ters as in Fig. ^. 



ppol 



E' H„ 




(b) log 





ft)K>gi„|H/ 



spol {b)log^JH|= 




H I E„ 



LJ 



(d) bg.JH/ 



FIG. 10: Near field intensity images in p polarization [(a) 
electric, (b) magnetic, the latter split into (c) x, horizontal 
component, and (d) z, vertical component, all of them in a 
logjg scale] in an area of 0.5x0.5 pm^ close to the fractal 
surface in Fig. ^ with 3 — 257 nm and £,l = 12.85 nm, where 
a strong localized optical mode is observed. Other parameters 
as in Fig. H. 



